Transverse Energy Density Fluctuations in the Color Glass Condensate Model 
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We calculate the transverse correlation of fluctuations of the deposited energy density in nuclear 
collisions in the framework of the Gaussian color glass condensate model. 



I. EVENT-BY-EVENT FLUCTUATIONS 

The event-by-cvcnt fluctuations of the transverse emis- 
sion pattern of hadrons in high-energy collisions of iden- 
tical heavy nuclei have recently attracted much interest 
experimentally [IHZ] and theoretically [51422). When av- 
eraged over collision events, the azimuthal angular dis- 
tribution of emitted hadrons around the beam axis is 
symmetric with respect to the plane perpendicular to 
the impact parameter vector b between the two nuclei 
[Au+Au or Cu+Cu at the Relativistic Heavy Ion Col- 
lider (RHIC) or Pb+Pb at the Large Hadron Collider 
(LHC)]. The event averaged angular distribution 

^ = ^( 1 + E-«(^) cos M' « 

where <j) p is the angle between and b, is therefore 
completely characterized by the even Fourier coefficients 
v n . The dominant coefficient, V2, is called elliptic flow. 

Owing to quantum fluctuations in the density distri- 
butions of the colliding nuclei and finite particle number 
effects on the distribution of emitted particles, the left- 
right symmetry is broken in individual collision events. 
The angular distribution can then be written in the form 

dN dN / ^ , , , , \ 
dV=^|l 1 + L^) C0S (^ + ^)J' ( 2 ) 

1 n—l 

where ip n describes the tilt angle of the "event plane" 
for each Fourier coefficient with respect to the reaction 
plane defined by the vector b. For even n, ip n is peaked 
around zero; for odd n, ip n is randomly distributed. The 
dominant odd coefficient, V3, is known as triangular flow. 
The event averages of the amplitude of the coefficients v n 
are found to be constant over a rather large pseudorapid- 
ity range (|7y| < 2) in Pb+Pb collisions at the LHC [5], 
indicating an approximately boost invariant origin. 

The main dynamical source of event-by-event fluctua- 
tions in the coefficients v n are believed to be nearly boost 
invariant fluctuations in the transverse distribution of 
the energy density at the beginning of the hydrodynamic 
expansion of the quark-gluon plasma formed in the nu- 
clear collisions. The geometric anisotropy of these fluc- 
tuations then translates in an anisotropic transverse col- 
lective flow pattern, which manifests itself in anisotropic 



particle emission. In the color glass condensate model of 
energy deposition there are two obvious sources of fluc- 
tuations in the deposited energy density. One is geomet- 
ric fluctuations of the position of nucleons in the collid- 
ing nuclei at the moment of impact, leading to trans- 
verse fluctuations in the density of field generating color 
charges. This mechanism has been studied widely and 
is usually described geometrically by the Monte-Carlo 
Glauber model [23] . The transverse correlation length of 
the fluctuations generated by this mechanism will be of 
the order of the nucleon radius. 

The other source of energy density fluctuations are 
fluctuations in the color field strength for a given density 
of color charges. This mechanism has not been investi- 
gated quantitatively up to now. The transverse correla- 
tion length generated by color field fluctuations will be 
dictated by the single scale governing the physics of the 
color glass condensate, the saturation scale Q s . Since 
Qj 1 is much smaller than the nucleon radius, the color 
field fluctuations can be expected to govern the micro- 
scopic structure of the transverse energy density fluctua- 
tions, which then is modulated on longer transverse scales 
by fluctuations in the nucleon density in the colliding nu- 
clei. 

Here we calculate the transverse correlation function of 
the deposited energy density in nuclear collisions in the 
framework of the Gaussian approximation to the color 
glass condensate originally proposed by McLerran and 
Venugopalan [2U [H] . In Section II we derive the two- 
point correlator of the energy density following the colli- 
sion of two color glass condensates. In Section III we eval- 
uate the resulting integrals and present numerical results 
for the correlation function for a representative choice of 
parameters. 



II. ENERGY DENSITY FLUCTUATIONS 

In the Gaussian random source approximation to the 
color glass condensate model of small-x gluon structure 
of atomic nuclei [HI [55] the probability distribution of 
color charge density p a (x) in the transverse plane is as- 
sumed to be of the form 

P[p a \ = cxp (-^2 / d 2 xp a ( X ) 2 ^j . (3) 
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Here /i 2 represents the area density of color charges in 
the colliding nuclei, and Q s — g\i is called the saturation 
scale, because it represents the scale at which the small- 
x evolution of the gluon density becomes nonlinear due 
to saturation effects [37J . Owing to the independent 
contributions of several nucleons to the color field, the 
Gaussian approximation is expected to provide a good 
description to the color source distribution in colliding 
nuclei at small x In the light-cone gauge, the Gaus- 
sian color charge distribution translates into a Gaussian 
distribution of transverse gauge field strengths. Here we 
will follow the work of Lappi [25] . 

To calculate the initial state density fluctuations 

( £ (x) £ (y)) - (e(x)) (e(y)) , (4) 

where x, y denote vectors in the transverse plane, we 
start from the expression for the deposited energy density 
of the gauge field given in Eq.(lO) in ref. [29] : 

e(x) = i*g(x)F£(x) + 2A" c (x)^' c (x) (5) 



with transverse vector indices i, j, m, n, ... = 1, 2. Follow- 
ing the collision, the field strength tensor in the region 
between the receding nuclei only receives contributions 
from the mixed terms, as the color field of each individ- 
ual nucleus is a pure gauge and the field strength tensor 
of each individual nucleus is thus zero outside the nuclear 
volume: 

i^.(x) = 5 / abc (^(l;x)^(2;x) 

+^(2;x)4(l;x)) (6) 
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A" c (x)^ c (x) = 9 -f abc fa>Vc (7) 
x A?(l;x)4(2;x)A«'(l;x)4(2;x). 

Here "1" and "2" denote the gauge fields carried by nu- 
cleus 1 and 2, respectively. The field correlator in the 
color glass condensate model is given by 



(4»(n;x)^( m ;y)) = i(yl?(n;x)4(m;y)) + \(A){m; y)A<}{n; x) 

d 2 p 



(27T) 



2 cos[p.(x-y)]^G(|p|) 



where G(|p|) is the Fourier transform of the function 



G(|x|) 



g 2 N\x\ 2 



1 — exp 



9 2 N i 
8tt 



V|x| 2 ln(A|x|) 



e(i-A|x|) 



with the IR cut-off parameter A. It is convenient to decompose the momentum quadrupole tensor as follows: 

P1+P2 r ,Px ~V\ 3 , 1 
PiPj = g ij + 2 ij +PlP2 CT y ' 

where a 1 , a 3 are the familiar Pauli matrices. We thus obtain 



(8) 



(9) 



(10) 



(^(n;x)A h j (m;y)^ = ^S mn 5 ab j (caa\p x (xx - y%)] cos\p 2 (x 2 - y 2 )] G(|p|) 
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+ cos[pi(a;i - yi)} cos[p 2 (x 2 - y 2 )} erf, G(|p|) 



■sintpi^x - yi)] sin[p 2 (a;2 - V2)] o\ 3 G(|p|) 



2 S mn S ab (SijD(x - y) + crfjE(x - y) - o-}jF(x - y)) = 5 mn S ab Sy(x - y). 



For later use, wc will note the values of the individual correlation functions D, E, F at the origin: 



d(o) = / (Si G (ipi) = ,^ G (i x i); E ^ = = °- 



(11) 



(12) 



The expression for D(0) diverges logarithmically for the function G(|x|) given in ([9]), if the gauge coupling g is taken 
as a constant. However, as pointed out by Kovchegov and Weigert [30) . the infrared divergence can be removed by 
including effects from the running of the coupling constant by means of the substitution 



ffV)5 a (l/|xf) 



(13) 
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in the exponent of The specific structure of this substitution, sometimes called the "triumvirate" structure of the 
running coupling, is motivated by the form of next-to-leading order corrections to the small-a; evolution of the BFKL 
kernel in the color dipole approach to parton saturation [3"T] . 

We first evaluate the expectation value of the deposited energy density: 



(e(x)) = ^/ obc / a ,^^(l;x)^(2;x),4? (l;x)A* (2; x) + A?(l; x)^(2; x)A? (2;x)A$ (l;x) 
+^fabcfa>b>cU?(l; x)^(2; x)A?' (1; x)4 (2; x) 



2 
2 

= yiV(iV 2 - l)D 2 (0), 

recovering the result given in Eq. (14) of ref. [29] , 

Next we evaluate the two-point correlator of the energy density: 



(14) 



<e(x)e(y)> 
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fabcfa>b>cfefdfe>f> d ({AUl;x)A b J2;x)A? (1; xU b - (2; x) 



+A Q (1; x)^ b (2; x)A?'(2; x)4 (1; x) + ^(1; x)A fc (2; x)A/(l; x)A$' (2; x)) 
x (> n (l; y)^{(2; j,)<(1; 2/ )<(2; y) + A e m (l; y)A{{2; t/)<(2; y)^'(l; y) 
+J 4 r e „(l; i / )4,(2; y)A( (2; y)" 



(15) 

We again make use of the fact that only correlators among fields in the same nucleus are non-zero, which allows us 
to suppress the labels 1 and 2: 

<e(x)E(y)) = ^fabofa'b'cfefdfe'f'd ((^° (x)A ? ' (x)^(y) <(y))(4(x)4 (x) ^ (y) < (y) 



+^"(x)A Q 

+(A?(x)^ 
+(A?(x)^ 



(x)^„ (y)^/' (y) ) (A) (x)4 (x)^ (y (y 
(x)4 e „(y)^' (y))(4(x)4 (x)A4(y)^' (y) 
(x)^ n (y)<(y))(^(x)^'(x)A/(y)^'(y) 
(x)A e m (y)A( (y) ) (a§ (x)j4"' (x)^ (y)< (y) 
(x)^( y )^'(y))(4(x)Af(x)^( y )A{'( y ) 
(x)A r e „(y)^;(y))( J 4 b (x)4(x)A{(y)A/;(y) 
(x)^(y)^' (y) ) (A b (x)4 (x) A* (y)< (y) 



+ (^(x)Af (x)^(y)^' (y))(^(x)4 (x)^(y)A£ (y))) 



5_ 
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fabcfa'b'cfefdfe'f'd ^ M a . 



(16) 



In the spirit of the Gaussian approximation we now factorize the correlators of four gauge fields into products of 
correlators among two gauge fields [35], e. g.: 

A^A? (x)A{(y)Ai(y)) = (^(x)^' «) (^(y)<(y)) + (^(x)^(y))(^'(x)A e m (y) 

4(x)A^(y))(^(y)^'(x)} (17) 

To proceed further we use the symmetry with respect to the color indices e' and /' to combine, e. g., the second 
and third term in the large brackets of (|T6|): 



Mo + M- 



3 = (A ? (x)A ? '(x)^(y)A{'(y))(4(x)4(x)(A/(y)A^(y) - A4(y)^'(y) 



(18) 



The second factor is easily shown to vanish: 

= (^(x)^(y))(4(x)<(y))(4(x)A/(y))(^(x)<(y)) - (m n) 

= $,/<W S jn (x.- y)S jm (x - y) + 5&//5& e ' - y)S,- m (x- y) - (m n) = (19) 

The same holds true for the forth and seventh term, M4 and M7. After considerable algebra, the fifth, sixth, eighth 



M 5 + M 6 + 



and ninth terms combine to 

M 8 + M 9 = ^/^/„^ c /e/d/e'/'«l(^(x)i4/(x)^(y)^;'(y) 

x( [^(x)4(x) - A$(x)4'(x)] [^(y)^(y) - ^(y)^'(y) 

~^fabcfa'b'cfefdfe'f'd (j^ae^a' f ^im^jn ~t~ ^af'^a'e^in^jm^j 
'• (dbffib'e''2[SinSj m — Si m Sj n ] + Sbe' Sf f2[Si m S j n — 



= |!^(^ - 1) [j>( x - y) » + £(x - y) 2 + F(x - y) 2 ] 2 , 



where we made use of the relation 
N 

fabcfa'b'cfa'bd — ~^~fab'd 



which follows from the Jacobi identity. 
(21) Finally we evaluate the first term: 



(20) 



M X = 9 -f a bcfa'b'cfefdfe'f>d ( A°(x)^' (x)^(y)^(y)} (^(x)^' (x)^(y)^' (y) 

= 9 -^N 2 (N 2 - 1) 2 D{0) 4 

+ 9 1 N 2 {N 2 - 1)D(0) 2 [D(x - y) 2 + £(x - y) 2 + F(x - y) 2 ] 

+ ^N 2 (N 2 - 1) [D(x - y) 2 + £(x - y) 2 + F(x - y) 2 ] 2 . 



(22) 



Combining these equations we finally obtain 

7g 4 



( £ (x) £ (y)) - (e(x))( £ ( y )) = y -N 2 (N 2 - 1)D(0) 2 [£>(x - y) 2 + £(x - y) 2 + F(x - y) 2 ] 
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iV 2 (7V 2 - 1) [£>(x - y) 2 + £(x - y) 2 + F(x - y) 2 



(23) 



Since the functions D,E,F always appear in the same combination, it makes sense to introduce the abbreviation 

K(x - y) = £>(x - y) 2 + £(x - y) 2 + F(x - y) 2 , (24) 
In terms of which the average deposited energy density and its fluctuation can be expressed as: 



so = {e) = 9 -N{N 2 -l)K{0)- 
Ae(x-y) 2 = ( £ (xMy))-(e(x))(e(y)) = ^ 2 (iV 2 -l) 



K(0)K( X -y) + -K(x-y) 2 



(25) 
(26) 
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III. EVALUATION OF INTEGRALS 



Next we simplify the integrals £>(x — y), E(x — y) and F(x — y). We abbreviate z = x y and z = |z|. We begin 
with D(z). 



D(z) 



We substitute 



^— % G(p) [ d<p cosfpzi cos cf>) cos(pz2 sin < 
o 47r Jo 



pdp 
4^2 



G(p) J d(f> ^ cos[p(zi cos + Z2 sin </>)] + cos[p(zi cos <fi — Z2 sin cj))] 
in the last term and introduce the notation z\ = zcosip, Z2 = zsini/;: 



G(p) / d0 cos[pzcos(</) — ip)] 



pdp 

~2tt~ 



G(p) J Q (pz) 



cPp 
4tt 2 



G(p) e*" = G(z), 



where we have used Eq. (3.715.18) from [53] . Similarly we obtain 



pdp f 2lT 

G(p) I d(f> cos(pz 1 cos <p) cos(pz2 sin </>) (cos 2 ^ — sin 2 . 
Jo 



dp 

10 4^ 2 j 
Using the same substitutions we find: 



G{p) j dcf> ^cos[p(zi cos 4> + z 2 sin (j>)] + cos[p(zi cos 4> — ?2 sin 4>)]^j cos(20). 



= — cos(2i/j) 



pdp f 27T 

G{p) I d(f> cos[pz cos(0 — -0)] cos(20) 



o 



pdp 



G(p)J 2 (pz) 



fo 2tt 

where we used Eqs. (3.715.18) and (3.715.7) from [52]. Finally, a similar calculation yields: 



F(z) 



pdp 



G{p) I d<j) sin(pz! cos 4>) sin(pz 2 sin (f>) cos <p sin <p 
o 47r Jo 



sin(2?/;) 



~2T 



G(p)J 2 ( P z). 



r 



(27) 



(28) 



(29) 



(30) 



(31) 



We conclude that the function K (z) only depends on the We have to evaluate integrals of the type 
distance z between the points x and y. 

We can express D(z) and E(z) 2 + F(z) 2 in terms of 
the following integrals: 



B n (x,z) = 



pdp 



2tt 



J Q (px)J n (pz). 



C n {z) 



pdp 



G(p)J„(pz), 



for n = 0, 2, namely: 
D(z) = C (z); 

E(z) 2 + F(z) 2 = C 2 (z) 2 . 
We rewrite the integrals as follows: 

C n (z) = J J d 2 xe-^G( X )J n (pz) 



(32) 

(33) 
(34) 



(36) 



(35) 



This integral can be evaluated for n = using formula 
(6.633.2) in [52]: 



f°° _ 2 2 e -(x 2 +Z 2 )/ic 2 f XZ\ 

J pdpe c p J (px)J (pz) = — 2 1 { - - i 



V2c 2 y 

(37) 



We are interested in the limit c — > 0, which means that 
we can apply the limit of Io(z) for large arguments: 



Io(z) -> 



'2ttz 



(38) 



(> 



This yields 



Bq(x,z) — lim -j— j= cxj) 



c->o Aire \J~kxz 
1 



(x - zf 
4c 2 



2?rz 



(5(x — z). 



(39) 



For n = 2 we use the recursion relation for Bessel func- 
tions: 



■h(z) = -Ji(z) - J (z), 
z 

and apply formula (6.512.3) from 



1 

dp Jo(px)J\(pz) — -6(z — x) 



with the convention 9(0) = 1/2. This implies: 



B 2 (x, z) = — ^8(z - x) 

■KZ Z 



2ttz 



5(x — z). 



(40) 



(41) 



(42) 



When we insert these results into the desired integrals, 
we find: 



C {z) = G(z); 



C 2 (z) = 



xdxG(x)-G(z). 



IV. NUMERICAL RESULTS 



(43) 
(44) 
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FIG. 2: The function Ci(z) for the selected parameters. 



Note that the result is independent of the value of A. 
The functions Cq(z) and C*2(z) are shown in Figs. [I] and 
[2] for these parameter values. 

For these parameters, the initial value of the deposited 
energy density is £o ~ 240 GeV/fm 3 . This very large 
energy density quickly decreases due to the longitudinal 
expansion and reaches much smaller values by the time of 
thermalization. What matters for us is not the absolute 
value of the initial energy density, but the relative size 
and spatial correlation of its fluctuations, Ae(x — y)/eo- 
This function is shown in Fig. [3] for the parameters listed 
above. 



O 3 
U 




0.0 



FIG. 1: The function Cq(z) for the selected parameters. 

We now evaluate the average energy density and its 
fluctuations for a choice of the parameters that is mo- 
tivated by the initial conditions at which thermal QCD 
matter is formed in heavy ion collisions at RHIC and 
LHC: 



Ql 
5 2 (M 2 ) 

Hl/x 2 ) 



( ff V) 2 = 2 GeV 2 ; 
3.785; 

16tt 2 
91n(l/(A 2 a: 2 )V 




(45) 



FIG. 3: The function Ae(z)/eo for the selected parameters. 



As the figure shows, the fluctuations of the initial en- 
ergy density are locally of similar magnitude as the en- 
ergy density itself and fall over distances of the inverse 
saturation scale, here assumed as Q^ 1 « 0.14 fm. This 
result is in accord with the intuitive picture of the field 
configuration immediately after the collision in the color 
glass condensate model, as a bundle of longitudinally 
stretching random color flux tubes with characteristic 
transverse width l/Q s - 
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V. SUMMARY 

We have calculated the initial energy density fluctua- 
tions in high-energy heavy ion collisions within the Gaus- 
sian color glass condensate model. These turn out to be 
very large with a transverse profile determined by the sat- 
uration scale Q s . A finite result is only obtained when 
the "triumvirate" running coupling is used, giving addi- 
tional support to the correctness of Eq. (13). The fluctu- 
ation probabilities thus derived can serve as input for any 
calculation aiming at the investigation of early fluctua- 
tions, in particular for calculations which study the fate 
of such fluctuations during thermalization. For example, 



it is possible to investigate the problem of event-by-event 
fluctuations in heavy ion collisions within the AdS / CFT 
paradigm using methods similar to those employed in 
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